Turbulence in a free surface 
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We report an experimental and numerical study of turbulent fluid motion in a free surface. 
The flow is realized experimentally on the surface of a tank filled with water stirred by a vertically 
oscillating grid positioned well below the surface. The effect of surface waves appears to be negligible 
so that the flow can numerically be realized with a flat surface and stress-free boundary conditions. 
The surface flow is unconventional in that it is not incompressible and neither energy nor enstrophy 
are conserved. Nevertheless, according to both experiment and numerical simulation, the second 
order structure function 5*2 (-R) scales essentially as for a three-dimensional system. However, the 
surface flow seems to be more intermittent. 
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Statistically stationary turbulence is readily produced 
in a tank of water by vertically oscillating a grid below the 
water surface. The properties of this three-dimensional 
volume turbulence have been studied extensively j|] . Par- 
ticles floating on the free surface are driven by the hor- 
izontal component of the turbulent flow below. While 
this is a classic problem, generally studied in an oceano- 
graphic context where the "particles" are micro organ- 
isms, such as phytoplankton or Lagrangian drifters 
[^,|| , controlled laboratory measurements of the relative 
velocity of floating particles do not appear to exist, de- 
spite a very large number of oscillating grid experiments 
which probe the fluid flow in the bulk and even near the 
surface . Such surface flows also seem to have received 
little attention by those interested in the fundamental as- 
pects of two or three dimensional turbulence [^|J|]. Our 
experimental and theoretical results demonstrate that 
these two-dimensional surface flows have unusual proper- 
ties that make them interesting in connection with scalar 
transport. The work opens the possibility of studying cer- 
tain aspects of passive scalar dynamics in divergent flows 
beyond the ones used for the Kraichnan model which are 
Gaussian and white in time Q . Their properties are also 
needed for a hydrodynamic modelling of the distribution 
of flotsam driven by underlying turbulence or chaos [^j. 
Our study is also related to recent work on quasi two- 
dimensional turbulence, which is probed by tracking the 
motion of floating particles The surface flows being 
comprised of particles constrained to the surface in both 
the experimental and numerical work, arc qualitatively 
different from the motion of the driving fluid. 

Although the fluid is incompressible, the inability of 
the floaters to enter the bulk flow assures that they will 
have a non-zero two-dimensional divergence. The com- 
pressibility of the floating particles causes them to clump, 
an affect seen in the numerical and experimental real- 
izations. The main differences between their motion and 
that of the underlying fluid comes from the possibility 
that the floaters can exchange kinetic energy and enstro- 
phy with the water molecules below. Hence, the energy 
and enstrophy of the floaters still would not be constants 



of the motion, even if the underlying fluid were invisid 
and undriven . With the above conservation laws ab- 
sent, one lacks the traditional dimensional arguments to 
estimate the scaling forms of the velocity and vortic- 
ity structure functions, 5*2 (i?) and S^{R). The structure 
functions are the main focus of our analysis and are de- 
fined as 



and 



S2{R) = ([(v(x + R) - v(x)) • R/Rf 



= ([w(x + R)-w(x)]2) 



(1) 



(2) 



where R is the separation between two points and both 
X — {x, y} and R are in a horizontal plane. The vortic- 
ity, of magnitude Lu{x,y), is perpendicular to that plane. 
A third quantity characterizing the surface flow is the 
dimensionless compressiblity coefficient C, defined as 



C=((V-vf)/((Vvf) 



(3) 



which lies between and 1 if the turbulence is isotropic 
0. As usual, (• • ■) denotes an ensemble average. 

The measurements were made in a square plexiglass 
tank filled with water to a height that was varied. The 
bars of the grid were PVC and were square in cross sec- 
tion with dimensions of 1 cm. The grid is also square and 
conforms closely with the square shape of the plexiglass 
box. To minimize flexing of the grid when it is oscillating, 
a number of the PVC bars are replaced by brass bars. For 
some of the measurements, a vertical post was placed at 
the center of the square grid to further suppress flexing, 
thereby reducing the amplitude of the surface waves. The 
apertures of the grid are also square, their dimensions 
being L—3.66 cm. The grid is supported by four vertical 
rods which run through pillow blocks, assuring its smooth 
vertical motion. All the parameters used in our set-up 
were typical of prior oscillating grid experiments ex- 
cept that the lateral dimensions of the tank were some- 
what larger. The grid is driven sinusoidally from above by 
a 1/4 hp motor. The vertical amplitude of grid motion 
was A = 1.25 cm. Whereas these parameters were varied 
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by only a small amount, the distance Z between the grid 
and the water surface (at zo=0) was rather widely varied. 
Reducing Z increases the Reynolds number of the turbu- 
lence on the surface and in the bulk close to the surface 
(zo > 0)- For all the measurements presented here, the 
oscillation frequency / of the grid was 4.5 Hz ± 0.3 Hz. 

The flow was measured using the technique of par- 
ticle imaging velocimetry (PIV). There one illuminates 
the seed particles with a (horizontal) sheet of laser light 
and tracks their motion photographically. A CCD camera 
was centered at a distance that was 23 cm from the near- 
est wall. It was mounted above the tank where we could 
capture images of the surface particles or those in the in- 
terior. The particle illumination was furnished by a Nd- 
Yag laser, with the beam spread into a horizontal sheet 
by a cylindrical lens. This sheet of light was typically 10 
cm wide and 0.5 cm thick. The floaters were particles 
of various diameters ranging from 10 to 50 [ivn. The in- 
terior measurements were made using neutrally buoyant 
polystyrene spheres of diameter 10 /zm. The images were 
acquired using a commercial particle imaging velocimetry 
apparatus [|ll| and the particle tracking was done using 
in-house software. To form a statistically reliable ensem- 
ble average, each measurement of S2{R) lasted several 
minutes and included ~ 200 image pairs. Each S2{R) 
was determined from roughly a million particle pairs. 
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FIG. 1. Log- log plot of the second moment of the (lon- 
gitudinal) velocity difference at the indicated distances 
(Z) between the surface and the oscillating grid. The os- 
cillating frequency of the grid and the surface-to-grid dis- 
tances, Z, are indicated. Other parameters are Z — 8.2 
cm: i?eA=80, lo=3.6 cm, Vrms = ^-27 cm/s; Z= 18.2 cm: 
i?eA=50, lo=l-5 cm, Vrms=0-7 cm/s. For both sets of 
measurements C=0.5 ± 0.1. 

The surface and bulk turbulence are characterized by 
the outer scale of the turbulence Iq and the Taylor 
microscale A. The outer scale is defined as the inte- 
gral of the normalized velocity autocorrelation function 
C{R) = {v{x + R) ■ t;(a;))/v^„j5, where the brackets in- 



dicate an average over points x,y and over many im- 
ages. Here f^m^ = {{vi{x,y) — {vi{x,y)))'^) where i is one 
of the two arbitrarily chosen directions in the horizon- 
tal plane. Separate experiments established that the tur- 
bulence is almost isotropic in the horizontal plane (see 
below). The Taylor microscale is defined here as A = 
V {''^rms) / {id-v{x) / dx)'^) . The Taylor microscale Reynolds 
number is Re\ = Vrms^/v, where v is the kinematic vis- 
cosity (in water i/— 0.01 cm^/s). The surface motion is 
further characterized by the vertical component of the 
vorticity w and by the dimensionless compressibility co- 
efficient C [see Eq. (|)]. 

Figure 1 is a log- log plot of 52 (i?) in the surface mea- 
sured for Z= 18.2 cm (closed circles) and for Z— 8.2 cm 
(open circles). The grid oscillation frequency for the two 
sets of measurements are indicated. The outer scale of the 
turbulence Iq at Z=18.2 cm and 8.2 cm are 1.5 and 3.6 
cm, respectively. The Taylor microscale Reynolds num- 
bers Rex are 80 {Z = 8.2 cm) and 50 [Z= 18.2 cm). 
Technical limitations blocked measurement of 5*2 at R 
less than 1 mm. The measurements in Fig. 1 suggest that 
52 (i?) is approaching its saturation value 2w^„^ at large 
R. The more interesting case is 52 at smaller i?, where 
one finds that S2{R) oc i?^, with ^=1.6 and 1.4 at Z= 
8.2 and 18.2 cm, respectively. The dissipative range at 
even smaller i?, where 52 (i?) — > i?^, was not accessi- 
ble in the laboratory experiment For these, and the 
other surface measurements, the dimensionless compress- 
ibility C defined above, was 0.5 ± 0.1. Reynolds num- 
bers in the experiment and the simulation are too small 
to develop a significant interval with algebraic scaling 
in 52 (i?) between the dissipative, 52 (i?) oc i?^, and the 
saturation range R » Iq, where S^iR) oc There- 
fore the extended self similarity (ESS) representation 
is used and a local scaling exponent D2.3{R) = 
dlogg 52(i?)/dlogg G3{R) for all R in the interval 0.1 cm 
< R < 7 cm is found to give 0.70 ± 0.10 (see upper 
panel of Fig. 2). For the ESS plots, the third order struc- 
ture function G^iR) = (|(v(x -f R) - v(x)) • R/i?|^) is 
needed. 

Since the surface is entirely driven by the grid- 
generated turbulence in the water, it is important to 
measure 52 (i?) below the surface. We have made such 
measurements at distances zq in the range 0.5 to 2 cm 
and for various grid-to-surface distances Z . At Z 9.6 cm, 
zo= 2.0 cm, the second moment is very well fitted to a 
power law form in the interval 0.1 cm < i? < 4 cm. In 
this interval C=0.9 ± 0.2, which is somewhat larger than 
the Kolmogorov value of 2/3 Other relevant param- 
eters for this bulk measurement are Re\— 100 and Iq— 
3.1 cm. The corresponding parameters measured on the 
surface are Re\= 40 and ^o= 3.5 cm. At all values of 
i?, the bulk value of the second moment is roughly an 
order of magnitude larger than 52 (i?) measured at the 
surface. A dimensional argument suggests the crossover 
between bulk and surface behavior should take place at 
a depth 5z below the surface given by 5z ~ v/v where 
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V is the bulk rms value of Vz- This gives 5z ~ 0.1 mm. 
The numerical studies indicate that the prefactor in 
this dimensional estimate can be as large as 10, making 
Sz ^ 1 mm, which is still too close to the surface for us 
to measure. 
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FIG. 2. Local scaling exponent, D2,3(i?), of an ex- 
tended self similarity plot of second order moment vs. 
absolute third order moment. Upper panel: The param- 
eters for this experiment are the same as in Fig. 1 at 
Z=18.2 cm. Lower panel: Comparison of D2,z{R) in the 
bulk at zo=0.5 (triangles) with that on the two surfaces 
at zo=0 and 1 (diamonds). Data were taken from a nu- 
merical simulation at Re\ ~ 100. Dashed lines in both 
panels are drawn at £'2.3(i?)=0.73 and 0.71, respectively. 



S^{R)R^ / S2{R)- At R equal to the outer scale of the 
turbulence, Zq — 3 cm, Q ~ 1, which suggests that it is 
probably too large to satisfy the zero-vorticity require- 
ment of wave turbulence theory, even though this theory 
successfully accounts for the motion of particles floating 
on a container of vertically oscillating fluid where capil- 
lary waves are excited [|l^,0 . 

The experimental observations are, however, in good 
agreement with a model for turbulence in a flat surface 
bounding a three-dimensional volume with fully devel- 
oped turbulence governed by the incompressible Navier- 
Stokes equation |]l3| . This fluid is in a computational 
box with periodic boundary conditions applied to the 
lateral edges but with the fluid particles on the upper 
and lower surfaces obeying free-slip boundary conditions, 
with Vz — dzVx = dzVy = 0. The equations were solved 
by means of a pseudospectral method with resolution 
Nx X Ny X Nz of 256 x 256 x 65. A statistically steady 
state of the turbulent flow was achieved by forcing in a 
wavenumber range comparable to the inverse of the box 
height (for more details see [0). Here Re\ ~ 100. 
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We consider the possibility that the surface velocity fluc- 
tuations originate from gravity and capillary waves gener- 
ated by the 4.4 Hz oscillation frequency of the grid. Even 
though the amplitude of the surface waves was rather 
small (A ~ 0.5 mm) it is comparable to the thickness 
of the boundary layer near the surface. The appropriate 
wave turbulence theory treats the nonlinear coupling of 
surface modes r2(fc) of different wavenumbers k (l4|. This 
theory is of Hamiltonian form, which requires that the 
(2D) velocity flcld in the theory is obtained from a veloc- 
ity potential, so that the vorticity of the surface particles 
is zero. Measurements of the (vertically-directed) vortic- 
ity of the surface particles show that it is quite large. 
Figure 4 is a log-log plot of the second moment of the 
surface vorticity fluctuations, 5'i^(i?), for a grid-to-surface 
distance Z= 8.5 cm. This function saturates at a rather 
small value of R and approaches the limiting value 2(w^). 
A dimensionless measure (Q) of the surface vorticity is 



FIG. 3. Deviations of the local extended self-similarity 
(ESS) slope from the classical Kolmogorov scaling for dif- 
ferent heights zq of the averaging plane and for different 
orders n = 2 to n = 6. All data points are calculated 
from 2 X 114 statistically independent samples at Re\ ~ 
100 and by averaging Dn.siR, zq) f^'" scales R between 
I877 and 4I77. 

The local ESS scaling exponent I?2, 3(.R) lies in the range 
between 0.67 and 0.74 in the interval 5 < R/t] < 100 
with the Kolmogorov length rj = (j/'^/e)^/^ (lower panel 
of Fig. 2) . This result is in good agreement with the 
experimental data also shown in Fig. 2 (upper panel). 
The simulation results show a plateau between I877 and 
41?7 for the local scaling exponents in both cases, at the 
surface and in the bulk. Intermittency corrections in the 
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surface seem to be stronger than in the bulk. In Fig. 3, 
a systematic investigation of these deviations from the 
classical Kolmogorov scaling with respect to zq, the ver- 
tical position of the horizontal plane of averaging is pre- 
sented. Our results indicate increasing deviations toward 
the free surface. The two-dimensional surface compress- 
ibility C in the calculation is 0.5, in excellent agreement 
to the laboratory-measured value of 0.5 ± 0.1. This value 
of C is expected for isotropic turbulence. 
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FIG. 4. Second moment of vorticity differences for the 
floating particles, at Z=9.6 cm. Here /=4.7 Hz, Re\=35, 
lo=1.7 cm, and t;f.ms=0-67 cm/s. Again C ~ 0.5. Inset: 
Model calculation of Si^{R) at the free-slip surfaces, with 
Rex ~ 100. Here C=0.5. 

The above model was also applied to the calculation of 
S^{R), and the results are presented in the inset of Fig. 4. 
Both the calculated and measured S^{R) are qualita- 
tively similar and show that this function rises rapidly 
with increasing R in these log-log plots and reaches its 
saturation value, 2(a;^) at a rather small value of R . 

In summary, we have described a laboratory experi- 
ment and a model calculation for particles that float on 
a sea of three-dimensional turbulence. Even though the 
surface particles move in a plane, there appears to be no 
basis for describing the floating particle system in terms 
of the laws of two-dimensional turbulence or of turbulent 
surface waves. Important characteristics of the floating 
particle system are its large compressibility, large vor- 
ticity and the absence of conservation laws for energy 
and enstrophy. The measurements and calculations pre- 
sented here indicate the existence of an inertial range 
where S2{R) scales as R^, with ( quite close to 2/3, as 
in 3D flows and in the inverse cascade^ of 2D ttirbulence. 
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